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The recent discovery of a large CP violating asymmetry in ^ 7r+7r~e^e~ mode has prompted us to seach for 
the associated Tr^Tr^e+e" decay mode in the KTeV-E799 experiment at Fermilab. In 2.7 x 10** Kl decays, 

one candidate event has been observed with an expected background of 0.3 event, resulting in an upper limit for the 
TT^TT^e+e" branching ratio of 6.6 x 10~^ at the 90% confidence level. 

PACS numbers: 13.25.Es, 13.30.Ce, 14.40.Aq 



The discovery in the KTeV experiment of the Kl 
7r+7r~e"'"e~ decay mode |^ later confirmed by KEK-E162 
and CERN-NA48 Hi together with the observation of a 
large CP violating asymmetry in the angular distribu- 
tion between the e+e~ and tt+tt" planes, has generated 
interest in searching for Kl — > 7r'^7r'*e+e^. This letter 
details the first attempt to detect this mode. 

The Kl tt^ ■k~ e~ decay can proceed at tree level 
via several amplitudes Q| . The main processes are a) the 
CP violating bremsstrahlung process in which the Kl 
decays into 7r"'"7r~ where one of the pions radiates a pho- 
ton which internally converts into an e^e~ pair, and b) 
the CP conserving direct decay of the Kl into a tt+tt^ 
pair plus an Ml photon. Other smaller amplitudes in- 
clude c) the CP violating direct decay of the Kl into a 
TT+TT^ plus an El photon and d) the "charge radius" am- 
plitude describing the Kl — > Ks transition by emission 
of a virtual photon followed by the CP conserving decay 
A'5 — > TT+TT^. This process is similar to the regeneration 
of Kl Ks by scattering from atomic electrons, both of 
which are proportional to the mean squared charge radius 



of the neutral kaon. In contrast, the Kl 
cay, because the pions are neutral, has no bremsstrahlung 
amplitude. In addition, the direct emission of a Ml pho- 
ton or an El photon is suppressed since gauge invariance 
and Bosc statistics require a tt'^tt'^ angular momentum of 
l>2§. 

A consequence of the suppression of the direct 
Ml and El emission processes and the absence of a 
bremsstrahlung process is that the Kl — > 7r'*7r°e^e" 
branching ratio is expected to be considerably smaller |^ 
than the Kl — > 7r^7r~e*"e~ branching ratio measured to 
be 3.2 X IQ-'^ §. The Kl ^ 7r°7r°e+e" is predicted to be 
dominated by the CP conserving "charge radius" ampli- 
tude with additional contributions from the CP conserv- 
ing E2 and CP violating M2 direct emission amplitudes. 
Theoretical estimates of the Kl Tr'^Tr^e+e^ branching 
ratio vary considerably from a vector meson dominance 
model estimate of 2.0 x 10"^ and chiral perturbation 
theory calculations of 1 x 10"*° g and 2.3 x 10"*" §]. 

An observation of Kl — > 7r'*7r"e"*'e~ would provide an 
indirect measurement of the Kl 7r'*7r'*7 branching ra- 
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tio, a mode which is difficult to detect directly because 
of the serious backgrounds ||Jl^- Observation of this 
mode would also allow investigation of chiral perturba- 
tion theory contributions to ©(p**) and 0{p^). The extra 
contributions of the full chiral perturbation calculations 
beyond that of the charge radius amplitude of Rcf. P] 
can complicate the extraction of this amplitude P,p|,|l 1| . 

A search for the Kl 7r°7r°e+e~ decay mode has 
been conducted by the KTcV collaboration using the 
1997 E799 data set. In the E799 spectrometer configura- 
tion, two almost parallel neutral beams passed through 
the fiducial decay region of the KTeV spectrometer 95 
m< z < 158 m downstream of a BeO production tar- 
get. The beams were produced by (3 — 5) x 10^^ 800 
GeV/c protons per minute delivered in 20 second spills 
incident at an angle of 4.8 mrad on a BeO target. The 
composition of the beams was mainly if^'s and neutrons 
with a small admixture of A's, and S's. An integrated 
total of 2.7 X 10^^ decays took place in the fiducial 
region while the trigger was live. The energies varied 
between 20 and 200 GeV with a mean of 70 GeV. 

The KTeV spectrometer components which were used 
for the analysis of these decays consisted of four stations 
of drift chambers, two upstream and two downstream of 
an analysis magnet which provided a pt kick of 0.205 
GeV/c for momentum determination of charged tracks. 
Downstream of the last drift chamber, in succession, were 
a charged particle trigger hodoscope which provided the 
first level trigger requirement of two charged tracks and 
a 3100 element pure Csl calorimeter for photon and elec- 
tron id and energy determination. Surrounding the decay 
region and the elements of the detector were veto coun- 
ters that detected charged and neutral particles that were 
not in the acceptance. The complete KTeV E799 spec- 
trometer configuration is described in Ref. [|l| 

The data were collected with a trigger which required 
two charged tracks in the drift chambers, four or more 
showers with energies > 1 GeV and total energy in the 
Csl calorimeter be greater than 30 GeV and no significant 
energy be present in any segment of the veto counters sur- 
rounding the decay volume and spectrometer elements. 

Events were selected off line that had six energy clus- 
ters in the Csl, two of which were associated with the 
e~^e~ pair. The decay vertices were determined by 
minimizing the of ^ fit which included the differences 
between the invariant mass of the secondary particles 
from the vr" decays (77 or e"'"e~7 in the case of modes 
with TT^'s undergoing Dalitz decay as discussed below) 
and TT^ mass together with the distance between the ver- 
tex and the tracks. The incorporated the measured 
uncertainties of the energies and positions of the photon 
clusters in the Csl and the uncertainties of the track 
position due to the drift chamber resolutions and mul- 
tiple scattering. The event reconstruction checked all 
possible pairings of the photons to determine the best tt" 
combinations. 



In addition, a large Kl tt^tt^tt^ sample (where tt^ 
indicates a Dalitz decay, 7r° e+e~7) was selected from 
events which had an e~^e~ pair and five Csl clusters unas- 
sociated with tracks. The Kl — > tt^tt^tt^ signal was com- 
pared to a Monte Carlo of this mode incorporating detec- 
tor acceptances, efficiencies and resolutions. Good agree- 
ment was observed in numerous kinematic distributions. 
Using the detector simulation of the Kl — > 7r°7r°e+e^ 
decays, a region of 0.493 GeV/c^ < M^o^o^+^- < 0.501 
GeV/c^ and {e+ e~ jjjj) < 1.5 x 10"'' GeVVc^ was 
determined for the 7r°7r°e+e~ search. The sig- 

nal region contained 86% of the Kl 7r"7r"e+e^ decays 
(using the Sehgal model for our generated signal events) . 
Here, (e+e~7777) was calculated relative to the vector 
from the center of the target to the decay vertex. 

The most serious background to the Kl — > 7r'^7r°e"'"e~ 
mode resulted from Kl Tr'^Tr'^Tr" decays in several ways: 
a) one photon from a tt" — > 77 decay converted inter- 
nally or externally into an e~^e~ pair while one of the 
other photons in the decay fails to produce a detectable 
energy deposit in the Csl or missed the Csl and the re- 
maining four photons formed a good tt'^tt*' pair; b) two of 
the six final state photons converted internally or exter- 
nally into an e+e^ pair and an and e~ were lost; c) 
one of the ir'^ decayed into an e"'"e~. In addition to the 
backgrounds from Kl 7r'^7r'^7r° decays, Kl — > tt^tt^ 
decays with an extra accidental Csl energy cluster con- 
tributed a background. Finally, other backgrounds due 
to Kl — > 7r''e'^e~7, — > tt^tt^j, and — > tt^tt'^jj 
decays were found to be negligible due to physics cuts 
discussed below and to small branching ratios. 

Special cuts were required to eliminate each type of 
background. The background due to Kl tt^tt^tt" de- 
cays in which two photons converted into e'^e~ pairs or 
one of the tt'^'s decayed into two e"'"e~ pairs was sup- 
pressed by rejecting events in which there were extra 
track segments upstream or downstream of the magnet. 
If a photon converted externally in the upstream spec- 
trometer material, the resulting e^e^ did not typically 
produce two reconstructed upstream track segments due 
to the small opening angle. If one member of the e"'"e~ 
pair had low momentum and was swept out of spectrom- 
eter acceptance by the analysis magnet, no extra track 
segments were produced downstream of the magnet ei- 
ther. In these cases, the invariant mass of the observed 
charged track and the other photon from the n'^ decay 
was close to the mass of the parent tt^ and the event was 
rejected by cutting on the 76^ mass. Specifically, the 
background from Kl — » tt^tt^tt'^ decays with two pho- 
tons from different tt^ converting to e~^e~ pairs and no 
extra track segments in the spectrometer was suppressed 
by rejecting events with two ej mass combinations with 
0.115 GeV/c^ < Me^ < 0.145 GeV/c^ and the two re- 
maining photons' mass within 3 MeV/c^ of the tt" mass. 

The background due to Kl — > tt^tt^tt" decays where 
one photon converted into an e^e~ pair and the other 
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photon missed the calorimeter were rejected by calculat- 
ing the missing photon energy and momentum assuming 
either 1) the missing photon went into the beam hole 
and had, in combination with the e~^e~ , the tt^ mass or 
2) the missing photon came from a Kl which had zero 

. If either, in the first case, M^+g-^^^^^ < 0.51 GeV/c^ 
and Pt (e"'"e~77777) < 0.0005 GeV^/c^ or, in the second 
case, the e^e^ plus missing gamma mass was greater 
than 0.065 GeV/c^, the event was rejected. 

To eliminate Kl — > Tr^Tr^Tr" decays where one photon 
converted into an e'^e" and the other photon was lost due 
to coalescing with another photon in the Csl, all observed 
clusters were tested to see if any cluster could be due to 
an overlap of two photons. The energy of each observed 
cluster in turn was divided in various proportions and the 
two resulting clusters were used in a fit of the event to 
the Kl — > 7r°7r°7r^ hypothesis. All possible combinations 
of all different clusters with different split fractions were 
combined with the other photons and the e"'"e~ pair were 
considered to see if any combination was consistent with 
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event. 



If the 7r°7r°7r^ invariant mass 



or the vertex x of a combination was consistent with a 
Kl — > 7r°7r°7r^ decay, the event was rejected. 

The background due to Kl tt^tt^tt" decays in 
which one ir^ decayed into an e~'"e~ pair (BR(7r° — > 
6+6")= 6.09 X 10~^ iQ) was eliminated by requiring 
that Me+e- < 0.10 GeV/c^. This cut also ehminated 
the background from a — > tt'^tt'^tt^ decay in which 
one of the tt^'s decayed into e'^e~e'^e~ and one of the 
e"^e~ pairs had low energy and went undetected. The 
background due to Kl "■^^'d decays with an acciden- 
tal cluster in the calorimeter was eliminated by requiring 
that Me+e-777 < 0.45 GeV/c^. All e+e~777 combina- 
tions were checked and the mass cut was applied to the 
one with the smallest for a '^'^t^d decay. 

After these cuts, the major remaining background, as 
determined from Monte Carlo simulations, was due to 
Kl — > tt'^tt'^tt^ decays. All other backgrounds were negli- 
gible. The contributions from all types of Kl tt^tt^tt^ 
backgrounds are shown in Fig. ^ as a function of pf and 
Mg+e-^^^^ for the equivalent of three E799 data sets af- 
ter all cuts. One event is observed in the signal region 
and the backgrounds outside the signal region are simi- 
lar in distribution and magnitude to the background ob- 
served in data when normalized to the E799 integrated 
flux. One event from the background Monte Carlo of 
three E799 data sets leads to an estimated background 
to the Kl TT^TT^TT^ signal mode of O.SS^q'ji events. 

Using the vertex fit, the mass of Mge-y^-y^ and of 
the e"'"e~7777 events from the triggers requiring two 
charged track and six cluster were calculated and the cuts 
mentioned above determined from Monte Carlo simula- 
tions of the signal and various backgrounds or from the 
Kl — ^ tt^tt'^tt^ data events were applied. Figure || shows 
the data remaining after applying the cuts. One event 
was observed in the signal box. Shown in Fig. ^ are the 



comparisons of the data and Monte Carlo simulations of 
the events vs. and Mg+g~-,^_,^. 
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FIG. 1 

Kl 7r^7r"e"^e~ pf vs. M^+^-^^^^ 
background events in the plot is equivalent to three times the 
actual integrated flux of the KTeV data. 
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FIG. 2. Events from the experimental data surviving all 
criteria for Kl -K^'n'^e'^ e~ decays vs. and M^+^-^^^^. 

An overall product of the acceptance and efficiency of 
0.25% for the signal mode was obtained from a Monte 
Carlo simulation of Kl Tr^Tr^e+e^ mode using the 
model of Ref. |^ with a chiral perturbation parameter 
Ws= 0.46 Since estimations of Ws from exisiting ex- 
perimental data are indirect and model dependent, we 
have studied the variation of our acceptance as a func- 
tion of Ws- We found that a simulation using a Ws= 0.9 
gives a minimum acceptance of 0.205% in the model of 
Ref. [|| (and a maximum upper limit for Kl ir^TT^e'^e^ 
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branching ratio). By comparison, a simulation using 
phase space for the decay gave an acceptance of 0.29%. 
We have used the lower acceptance 0.205% for calculat- 
ing the single event sensitivity using in determining a 
Kl TT^ Ti'^ e'^ e~ branching ratio upper limit. 

Kl ^ TT^TT^TT^ events with a missing 7 down the beam 
hole (which could be reconstructed as Kl tt^tt^tt^ 
using the technique for treating beam hole photons dis- 
cussed above) were used as a normalization mode since 
they had the same trigger and analysis cuts as the sig- 
nal mode. Figure ^ shows the mass distribution of these 
events compared to the Monte Carlo simulation. A dis- 
tinct Kl peak is reconstructed for the data with a res- 
olution matching the simulation of this type of events. 
The disagreement between the level of data and Monte 
Carlo of the normalization events has been folded into 
the systematic error in the single event sensitivity. 
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a)Mg+g~ and b) pt distributions of events surviving cuts. 
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FIG. 4. The mass distribution of Kl 
the photon from the n" — > e'^e~j decay passes into one of 
beam holes in the Csl calorimeter and is reconstructed as- 
suming that its trajectory passes through the center of the 
given beam hole. The solid dots are the data and the his- 
togram is the Monte Carlo simulation of these events. 
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A single event sensitivity of 1.71 ± Q.06{stat) ± 
0.19{syst) X 10~^ was calculated using the ratio of the 
Kl — > 7r"7r°e+e^ signal efficiency to the Tr^vr^vr^ 
normalization sample efficiency, the Kl — > tt'^tt^tt^ 
branching ratio of 7.59 x 10"'^ and the observed 
2234 Kl ^ tt^tt^tt^ normalization events of the type 
described in the previous paragraph. The observation of 
one event in the signal box resulted in an upper limit 
for BR{Kl "> Tr°Tr°e+e-) of 6.6 x 10"^ at the 90% con- 
fidence level. Using the ratio of the BR(i^i, 7r'^7r'^7) 
to BR{Kl 7r07r°e+e-) of «50 §, our result is a fac- 
tor of 10 more sensitive than the present upper limit of 
5.6 X 10~^ from a direct search [|l^ for Kl Tr'^n^'-f. 

In conclusion, we have made the first attempt to detect 
the Kl 7r°7r"e+e^ decay mode. We have established 
an upper limit of 6.6 x 10~^ at the 90% confidence level 
(using aws= 0.9 in the model of Ref. |^) which excludes 
the vector dominance model of Ref. but does not quite 
reach the level of expected branching ratios predicted by 
the more recent models of Refs. 
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